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ABSTRACT 


In  this  paper  we  consider  the  time  to  first  failure  of  a  parallel  system 
in  which  the  failure  and  repair  rates  of  components  depend  on  the  state  of 
the  other  components  as  well.  A  back-up  unit  with  a  random  life  time  is 
employed  whenever  all  the  components  of  the  system  are  down.  The  system 
fails  when  all  the  components  of  the  system  and  the  back-up  unit  are  down. 
The  first  moment,  the  Laplace  transform  and  the  probability  distribution  of 
the  time  to  first  failure  of  this  system  are  obtained.  Sufficient  conditions 
under  which  this  distribution  has  the  new  better  than  used  (NBU)  and  an 
exponential  limit  property  are  given.  Special  cases  with  phase  type  and 
deterministic  back-up  unit  lifetimes  are  also  considered.  These  results 
extend  the  results  of  Ross  and  Schechtman  (1979). 

Key  words:  Dependent  parallel  system,  first  failure  time,  back-up 
unit,  safety  periods,  new  better  than  used  (NBU)  distributions,  exponential 
limit  theorem,  phase  type  distributions. 


INTRODUCTION 


In  many  parallel  reliability  systems  with  maintained  components,  a 
system  failure  may  not  occur  immediately  upon  the  failure  of  all  of  its 
components;  rather  a  system  failure  may  occur  only  after  all  of  its  com¬ 
ponents  have  been  down  continuously  for  a  fixed  or  random  time  period. 
This  time  period  is  called  the  safety  period.  This  safety  period  may  be  the 
lifetime  of  a  back-up  unit.  For  example  consider  a  system  serviced  by  two 
AC  power  sources.  When  both  these  power  sources  fail,  a  DC  power  source 
(battery)  with  a  fixed  or  random  lifetime  is  used  to  supply  power  to  the 
system:  the  battery  here  serves  as  a  cold-stand-by  unit,  providing  a 
safety  period  for  the  AC  power  sources  to  be  repaired.  On  the  other 
hand,  the  safety  period  may  itself  represent  the  time  required  to  cause 
actual  damage  to  the  system  when  all  of  its  components  are  down.  For 
example  consider  a  vessel  accommodating  a  constant  heat  source  cooled  by  a 
system  of  two  redundant  blowers  (Dunbar  (1984)).  In  such  a  system  it  is 
usual  to  consider  a  system  failure  as  structural  damage  of  the  vessel  due  to 
excessive  temperature.  Suppose  the  normal  temperature  of  operation  be 
much  lower  than  the  maximum  that  the  vessel  has  been  designed  to  with¬ 
stand.  In  this  case  failure  of  both  blowers  is  not  sufficient  to  cause  im¬ 
mediate  failure  of  the  system.  The  system  will  fail  with  a  delay  in  time 
depending  on  the  magnitude  of  the  heat  input,  system  heat  capacity  and  the 
difference  between  the  operating  temperature  and  the  temperature  required 
to  cause  damage  to  the  vessel.  Furthermore,  the  system  will  fail  only  if 
neither  of  the  blowers  is  repaired  (or  replaced)  before  damage  occurs. 
Similar  examples  in  the  nuclear  and  perishable  food  industries  are  discussed 
in  Ross  and  Schechtman  (1979). 
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In  this  paper  we  consider  such  a  parallel  system  in  which  the  failure 
and  repair  rates  of  components  depend  on  the  state  of  the  other  components  as 
well.  A  back-up  unit  with  a  random  lifetime  is  employed  whenever  all  the 
components  of  the  system  are  down.  The  system  fails  when  all  the  components 
of  the  system  and  the  back-up  unit  are  down.  A  detailed  description  of  this 
model  is  given  in  Section  1.  The  first  moment,  the  Laplace  transform  and  the 
probability  distribution  of  the  time  to  first  failure  of  this  system  are  obtained 
in  the  same  section.  Section  1  also  contains  sufficient  conditions  under  which 
the  first  failure  time  has  the  new  better  than  used  and  an  exponential  limit 
property.  Special  cases  with  phase  type  and  deterministic  safety  periods  are 
considered  in  Sections  2  and  3  respectively. 


1.  DEPENDENT  PARALLEL  SYSTEMS  WITH  SAFETY  PERIODS:  THE  MODEL 


Consider  a  parallel  system  consisting  of  n  components.  Each  component 
alternate  between  intervals  in  which  they  are  up  (i.e.  working)  and  in  which 
they  are  down  (i.e.  failed).  Let  Z j ( t )  be  the  state  of  component  i  at  time 
t :  i .  e  Zj(t)  =  0  if  component  i  is  up  at  time  t  and  Zj(t)  =  1  otherwise.  Thus 
Z(t)  =  (Z^(t),  Z2U), . . .  ,Zn(t))  describes  the  performance  of  the  components 
at  time  t.  It  is  assumed  that  the  vector  performance  process  Z  ^{Z(t),  teR+} 
is  Markov  on  the  state  space  S  =  {0,  1}n,  and  if  {Zj,  ieN+}  is  the  embedded 
process  of  Z  embedded  at  its  transition  epochs,  1 1 Z j  -  Z j _ ^  1 1  =  1,  a.s.,  where 


for  a,  b  e  S,  lla-bll  s  1  (a-  -  b-)2.  That  is  at  any  transition  epoch  of  Z 

i  =  1  1  1 

either  one  repair  or  one  failure  of  a  component  takes  place.  We  call  such  a 
orocess  multivariate  binary  birth  and  death  process.  Note  that  if  all  n 
components  are  separately  maintained  and  the  failure  (respectively  repair) 
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rate  of  component  i  is  Aj  (respectively  (jj),  the  vector  performance  process 
Z  of  such  a  system  satisfies  the  above  condition  (e.g.  Chiang  and  Niu 
(1980),  Ross  and  Schechtman  (1979)). 

Our  system  is,  however,  dependent  in  the  sense  that  the  failure  and 
repair  rates  of  component  i  can  be  dependent  on  the  state  of  the  other 
components.  Suppose  the  vector  performance  process  Z  is  in  state  a  (re¬ 
spectively  b)  such  that  a ^  =  0  (respectively  bj  =  1).  The  failure  (re¬ 
spectively  repair)  rate  Aj(a)  (respectively  M , ( b ) )  of  component  i  may  depend 
on  a  (respectively  b)(e.g.  Ross  (1984),  Schechner  (1984)).  Such  a  generali¬ 
sation  is  needed  when  the  working  components  share  the  overall  workload  or 
when  the  repair  facility  has  limited  capacity.  Let  Q  be  the  infinitesimal 
generator  of  Z.  Two  examples  with  n  =  3  for  (i)  independently  maintained 
and  (ii)  dependent  system  with  equal  load  sharing  of  a  constant  load  and 
proportional  failure  rate  are  shown  in  Figures  1  and  2. 

When  Z  is  in  state  0  =  (0,0,0, ...  ,0)  (i.e.  all  n  components  are  working) 
we  say  that  the  system  is  in  perfect  condition.  The  system  is  said  to  be 
down  whenever  Z  takes  the  value  \  =  (1,1,..., 1)  (i.e  all  n  components  are 
failed)  and  it  is  said  to  be  u£  otherwise  (i.e.  at  least  one  component  is 
working).  Whenever  the  system  reaches  the  down  state  an  emergency 
back-up  unit  is  brought  in  which  provides  a  random  safety  period.  If  none 
of  the  n  failed  components  are  repaired  during  this  safety  period,  the 
system  fails  at  the  end  of  this  period.  Otherwise  at  the  time  of  the  first 
repair  after  the  system  became  down,  the  emergency  back-up  unit  is  removed 
and  the  system  is  placed  in  normal  operation. 

Let  0<R^<R2<.  •  .(respectively  0<S-1<S2<  •  •  • )  be  the  consecutive  time 
epochs  at  which  the  process  Z  enters  (respectively  exits)  state  2  starting 
from  some  initial  state  Z(0)  t  Then  one  sees  that  0<R.<S.<R2<S2<  •  • . . , 


a.  s 


Infinitesimal  Generator  Q  for  a  Separately 
Maintained  System  with  3  Components. 


(LVo)  U'o'i)  (o'l'l)  U'O'O) 


Figure  2:  Infinitesimal  Generator  Q  for  a  Dependent 

System  with  Equal  Load  Sharing,  Proportional 
Failure  Rate  and  Independent  Repairs. 


Let  XQ  =  Rr  Xj  =  Rj+-j-Sj  and  V.  =  Sj  -  Rj,  i  =  1,2,...  Note  thar  Yj  is  the 
time  required  for  the  first  repair  after  the  i-th  entrance  of  the  process  Z 
into  the  state  1_-  Thus  (Yj)*_^  is  a  sequence  of  independent  and  identically 
distributed  (i.i.d)  exponential  random  variables.  It  can  be  easily  verified 

Oo 

that  (Xj)j_Q  is  a  sequence  of  mutually  independent  phase  type  random  variables 
(e.g.  Neuts  (1981))  with  (X|)*=1  having  identical  distributions.  Now  let  (T.)°°=1 
be  the  sequence  of  the  length  of  the  safety  periods.  That  is,  at  the  entrance 
of  Z  into  state  2  for  its  i-th  time,  the  emergency  back-up  unit  will  provide  a 
safety  period  of  length  T..  Then  the  time  to  first  system  failure  VT 
is  given  by: 

VT  =  inf  {R.  +  T.:  R.  +  T.  <  S-,  i=1,2,....|  (1) 

In  the  remainder  of  this  paper  we  will  assume  that  (Tj)j_^  is  a  sequence  of 
independently  distributed  random  variables  that  are  independent  of  Z.  Note 
that  V-j.  is  nondecreasing  in  Rf  and  Tj.  Then  from  (1)  one  easily  sees 

Lemma  1:  For  X.'  S\  Xj',  i  =  0,1,2,...,  Y'  f' 1  Y !',  and 
V  |l  TV  ,  i  =  1,2,...,  one  has 

VT'  f  VT"  ' 

where  for  any  two  random  variables  W1  and  W",  W1^ 

(respectively  ^)  W"<=>P{W’  >  t}  S  (respectively  =) 

P{W"  >  t}  for  all  t. 

It  is  of  interest  to  obtain  the  mean  and  the  probability  distribution  of  Vj .  In 
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Section  1.1  we  will  first  obtain  these  results  with  no  specific  assumptions  regarding 
the  probability  distributions  of  (Xj)*_q  (Yj)*=1  and  (T. . 


1 . 1  General  Results  for  Single  Component  Systems. 

Consider  a  single  component  system  where  (Xj)*.-,  (Yj)*^  and  (T.)°°=1 
are  the  up,  down  and  safety  periods  respectively.  Then  the  time  to  the  first 
system  falure  V-p  as  defined  in  (1 )  can  be  rewritten  as: 


i-1 

V  =  inf  {  I  (X .  +  Y  )  +  T  :  T  <  Y. ,  i=1 ,2, . . . }  ,  (2) 

1  j=0  til 

where  Yq  =  0  w.p.1.  Let  Fq,  F,  G  and  H  be  the  probability  distribution 
functions  of  the  mutually  independent  random  variables  XQ,  (Xj)*^,  (Yj)*^ 
and  (T-( respectively.  Define  K+1  to  be  the  smallest  integer  for  which 
Tj  <  Yj  is  satisfied.  That  is 

K  +  1  =  min  {i:  Tj  <  Yj,  i  =  1 ,2,  .  .  . } .  (3) 

00  00 

Since  (Tj)j_^  and  ( Y . ) j _ ^  are  mutually  independent  sequences  of  i.i.d 
random  variables,  K  has  a  geometric  distribution  on  {0,1,2,....} 
with  paramete  P{Tj  <  Y{ } .  So 

P}K  =  i}  =  p(1-p)',  i  =  0,1,2 . 

and 


E { K }  =  (1  -  p)/p  , 
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where 

p  =  P{T.  <  V.}  =  J  (1  -  G(x))dH(x), 

1  1  0 

and  E  is  the  expectation  operator.  Now  combining  (2)  and  (3)  one 
obtains 


VT  =  X0  *  Yj*>  +  TiUl 


with  the  usual  convention  that  2  Xj  -  0  for  b  <  a  and  where 

a 


st  st 

Y?  =  Y„Y.<T.*nd  Y'Vj  ,t 


i  -  •  jlY,  <  T.  ~~  -j*  =  TilT.<  Y,'  i  = 


Note  that  for  j  ^  K,  Y*  =  Y.  and  TK+1  =  T*+r  The  probability  dis¬ 
tribution  functions  G*  and  H*  of  (Yj)*.,  and  (Tp*=1  respectively  are 
given  by 

i  x 

G*(x)  =  tt;  J  (1  -  H(y))dG(y)  ,  x  >  0 
I'P  o 

and 

H*(x)  =  ~  J  (1  -  G(y))dH(y) ,  x  >  0. 

p  0 


Taking  the  expectations  on  both  sides  of  (4)  one  obtains 


(4) 


E{Vt}  =  E{X0}  +  ^  { E  {X }  +  E{ Y*} }  +  E{T*}  , 


(5) 


where  E{X}  =  E{Xj},  i  =  1,2,. 

E{Y*}  =  E{Y*}  =  J  y(1  -  H(y))dG(y) , 

1  1  P  0 

and 

1  * 

E{T*}  =  E{Tf}  =  ^  J  y  (1  -  G(y))dH(y),  i=1,2, - 

1  p  0 

Furthermore,  taking  the  Laplace  transform  one  has 

L{Vy}  =  L{Xq}(  - ? - \  L{T*}  ,  (6) 

U  (  1-(1-p)L{X}L{Y*}  J 

where  for  any  nonnegative  random  variable  W,  L{W}  =  E{exp(-sW)},  Re(s)  >  0 
Hence 

P{VT  ^xl  =  Fn@  (  p  I  (i-p)jF(i)  ®  G*(j))  ®  H*(x),  x  >  0,  (7) 

T  0  j=0 

where  F^(x)  =  1,  x  >  0,  F^  is  the  j  (£1)  fold  convolution  of  F  with  itself 
and  ®  stands  for  convolution.  G*^  is  similarly  defined.  Equation  (6)  is  a 
generalisation  of  Equation  (2)  of  Ross  and  Schechtman  (1979).  When  Fq,  F, 

G  and  H  are  all  absolutely  continuous  with  support  (0,  »)  one  may  use  either 
the  Laguerre  transform  (Keilson  and  Nunn  (1979),  and  Sumita  (1981))  or 
the  generalised  phase  type  (Shanthikumar  (1985))  method  to  compute 
P{Vy  $  x}  given  in  (7).  However,  it  seems  appropriate  that  we  also  develop 
approximations  and  bounds  for  VT  that  are  easily  computable.  With  this  in 
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mind  we  will  first  provide  sufficient  conditions  under  which  Vy  has  the  new 
better  than  used  (NBU)  property.  A  nonnegative  random  variable  W  is 

st 

NBU  iff  P{W  >  t  +  s}  ^  P{W  >  t}  P{W  >  s}  or  equivalently  W  i  wiw>t-t  *  >  0. 

Theorem  2:  Assume  that  Xq  and  are  NBU,  (Vj)  °°=1  are  exponentially 

st 

distributed  and  Xq  £  Xj|x  >u"U'  u  >  *  =  1/2,  .  Then  Vy  is  NBU. 

i 


Proof:  In  order  to  establish  that  Vy  is  NBU  we  need  to  show  that 


st 


VT  -  VTIVT>t  't'  1  >  °* 


(8) 


Suppose  Vy  >  t.  Then  there  are  three  cases  one  can  consider.  They  are: 


Case  I 


0  <  t  <  R 


1 


In  this  case,  from  (4)  one  sees  that 


K 


VTICase  ,  -*  =  XOIX0>t-‘  (Xj  +  Vf>  *  TW 


(9) 


St 


Since  Xq)X  -t  S  Xq,  comparing  (4)  and  (9)  one  sees  that  (8)  is 
satisfied  in  this  case. 


Case  II:  R£  <  t  <  S£,  Vy  >  t  and  T£  >  t-R£  for  some  £(£  1). 


Let  u  =  t-R£  and  define 


st 

T '  =  T 


1  1  £IT„  >u'U'  Tj'  =  T£+j-1 '  j  =  2'3, 


'1 


sl+ j-r  1 


U  f  f  ,  •  .  •  f 


£IYjL>U  '  j  " 


and  X'  =  x£+j.1r  j  =  1/2, - -  and  set  =  0. 


Then 


i-1 

=  inf  {  1  (X*  +  Y.')  +  T.1  :  T.‘  <  Y.', 

T I  Case  II  j_q  J  J  i  i  i 


(10) 


i  =  1,2,-..}  . 


st 


Since  (T^ )j_-1  is  NBU  and  (Yj)*=1  is  exponential  one  sees  that  T*  ^  T.  and 


J  pi  'TP 

st  st 


Y.1  =  Y.,  j  =  1,2, _  Furthermore,  X1  =  X^,  j  =  0,1 _  Then  comparing  (2) 

and  (10),  from  Lemma  1  one  sees  that  (8)  is  satisfied  in  this  case  too. 


Case  III:  <  t  <  R£+1#  V-p  >  t  for  some  (£  £  1). 

Let  u  =  t-S£  and  define 
st 


*0  =  >u'U'  *j  X2  +  j'  1  ~  1,Z' - ' 

Yr  Vj-v  Tj'=  T£.j-v  i  =  -1'2 . and  Yo s  o. 


Then 


st 


i-1 


VT.Case  l..-t  =  inf^  j=IQ  +  Y/>  +  Ti  :  Ti  <  Yi  ' 


(11) 


i  =  1,2,...}. 


st 


st 


st 


Since  Xn  £  X0IV  -u  one  sees  that  X.  £  X1,  j  =  0,1,2,...,  Y.‘  =  Y.  and 

st  u  x  ^  u  J  i  J  J 

T.' =  T .,  j  =  1,2, _  Then  from  Lemma  1  and  Equations  (2)  and  (11)  one 

sees  that  (8)  is  satisfied  in  this  case  as  well. 
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Now  a  routine  unconditioning  with  respect  to  these  three  cases  (e.g.  Shanthikumar 
(1984))  completes  the  proof. 


The  above  theorem  appears  to  be  weaker  than  Proposition  1  of  Ross  and 
Schechtman  (1979).  As  we  will  see  later  in  Section  1.2,  application  of 
Theorem  2  to  the  dependent  parallel  system  results  in  a  theorem  that  is  stronger 
than  Proposition  2  of  Ross  and  Schechtman  (1979).  When  T.  is  degenerate 
(say  equal  to  a  constant  A),  one  can  use  the  proof  of  Proposition  1  of  Ross 
and  Schechtman  (1979)  and  the  above  proof  to  establish 

Theorem  3:  When  Xq  is  NBU,  T.  =  A  is  a  constant  and  for  all  u  >  0, 
st 

X0  ~  XilX.>u  *u'  '  *  VT  is  NBU‘ 

I 

Note  that  Theorem  3  is  a  slight  generalisation  of  Proposition  1  of  Ross  and 
Schechtman  (1979).  Theorems  2  and  3  can  be  used  to  bound  higher  moments 
of  V-j.  (e.g.  Proposition  3  and  Corollary  1  of  Ross  and  Schechtman  (1979)). 

Next  we  will  provide  an  exponential  limit  theorem  that  can  be  used  to  approximate 
the  distribution  of  V-j..  Consider  the  sum 

Ky 

U  =  I  (X.  +  vY) 

*  i=1  1 


where  is  a  geometric  random  variable  with  parameter  p^  and  support 

{0,1,2,....}, 


Pv  =  S  (1  •  G(yx))dH(x) 


and 


yX  =  Y. 


ilY.  <yT. 

l  i 

for  some  y  >  0.  Then 

O-P  J 


/  i  -  1/2 


E{uv}  = 


E{Xj }  +  ECYjIYj  <  v  T.} 


Suppose  P{T.  >  0}  =  1  and  Y.  has  support  (  0,  »)  and  is  absolutely 

continuous.  Then  as  y  ■*  °°  one  has  p^  -»  0,  E { U ^ }  -+  ®  and  L{yX)  -►  L{Y.}. 

Hence,  using  an  analysis  similar  to  that  of  Shanthikumar  and  Sumita  (1983). 

Theorem  1.A4,  one  obtains 

P{l>Y  i  t  E{Uy}}  -►  1  -  e'\  t  >  0  as  y  -»  °° 

Therefore 

Theorem  4:  When  G  is  absolutely  continuous  with  support  (0,  ®),  H(0)  =  0 
and  0  <  p  =  P{Yj  >  Tj}  <  e,  there  exists  an  M(e)  such  that  M(e)  -*•  0 
as  e  -»  0  and 

IP{VT  ^  x}  >  Fq  0  B  ®  H*  (x)l  <  M(e)  ,  x  >  0 
where 

B(x)  =  1  -  exp  {  - ^ -  }  ,  x  >  0. 

(1-p)(E{X}  +  E{ Y*} ) 


In  the  spirit  of  the  above  theorem  one  may  therefore  approximate  P{VT  i  x 


for  small  values  of  p  by 


P{VT  i  x}  =  FQ  ®  B  ©  H*  (x),  x  >  0 
Two  simpler  approximations  of  the  above  form  (for  small  values  of  p)  are: 

P{Vy  S  x}  =  BQ  ®  H*  (x),  x  >  0,  (12) 

and 

P{VT  i  x}  =  1  -  exp  {  x/E{Vt}};  x  >  0,  (13) 

where  EfV^}  is  as  given  in  (5)  and 

Bq(x)  =  1  -  exp  (x/(E{X0)  +  (^){E{X}  +  E{Y*} } )} ,  x  >  0. 

1 .2  Results  for  Dependent  Parallel  Systems 

00 

In  the  remainder  of  this  paper  we  will  assume  that  the  sequences  (X.)j=Q 

and  (Yj)°l-j  are  as  defined  for  the  dependent  parallel  system  governed  by 

the  vector  performance  process  Z  with  infinitesimal  generator  Q.  We  will 

00 

first  give  the  means  and  the  probability  distributions  of  Xq,  (Xj ) j_-|  ■  Partition 
the  Q  matrix  such  that 


S'  1 


1  hQ)  -m 


where  Qq  is  the  infintesimal  generator  of  the  lossy  process  of  Z  defined  over 
the  state  space  S'  ;  S  •  1  (Keilson  (1979)),  jjQ)  is  the  vector  of  rates  at 
which  the  process  Z  leaves  the  state  2-  The  rate  from  state  2  t0 
(1 ,1 , . . .  ,1 ,0,1 , . . .  ,1 )  when  the  0  is  in  the  i-th  position  is  MjQ)'  ■  =  1»2,...,n. 
All  other  rates  are  zero.  So 


. . -I--*.  -  1.  i«.*^  ’j.  V.  'J  -■-»  ‘j-  ’j  ^ 
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n 

M  =  I  M:(1)  =  H(1).  V. 
i=1  1 

For  the  example  corresponding  to  Figure  2,  jjO)  =  (0,0,0/0,n3,p2,p1 ), 

M  =  (p^  +  p2  +  M3)  and  QqV  =  (0/0,0/0/X;^/^)1  •  Then  (e.g.  Keilson  (1979), 
Neuts  (1981),  Shanthikumar  (1985)), 


E{X0}  =  aQ"1!1,  E{X}  =  b  Q’1  I1 


Fq(x)  =  a 


00 

2 

m=0 


Pme~AxfAx)m 

m! 


V,  x  >  0 


F(x)  =  b 


00 

1 

m=0 


Pm  e~Ax(Ax)m 
m ! 


V,  x  >  0 


where  a  is  the  initial  probability  distribution  of  Z(0)  and  b  is  the  probability 
distribution  of  Z  immediately  after  its  exit  from  state  That  is 
b  ( 1,1, ...  ,1,0,1, ...  ,1)  =  Pj(2)/M,  i=1,2,...,n  and  all  other  probabilities  are 
zero:  i.e.  b  =  jjO)/M.  Similarly  if  the  system  starts  from  perfection 
P{Z(0)  =  0}  =  a(0,0,...,0)  =  1.  Furthermore, 


A  =  max  {-Q0(i,i)} 

p  = 1  4  % 

and  I  is  an  identity  matrix  of  appropriate  dimension.  Also  one  can  easily  see 


that 
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E  { V  j }  =  1/M 

and 

G(x)  =  1  -  exp(-Mx),  x  >  0  . 

We  will  next  provide  sufficient  conditions  under  which  Theorem  2  can  be 
applied  to  this  case  of  multivariate  binary  birth  and  death  process  Z  defined 
on  the  state  space  S  =  {0,  1}n.  For  each  state  x  e  S  define 
W(x)  =  {i:  X.  =  0,  i  =  1,2,  ...,n}  and  W(x)  =  {1 ,2, . . .  ,n}-W(x) .  That  is  W(x)  is 
the  set  of  working  components  in  state  x.  Then  Shaked  and  Shanthikumar  (1985) 
has  shown  that  if 

IZ(0)  =  0  and  for  all  x,  y  e  S  such  that 
W(x)  c  W(y),  \j(x)  S  XjCy),  i  e  W(x) 
and  Mj(x)  £  >  i  eW(y), 

st 

then  Xq  is  NBU  and  Xq  t  >u'u<  u  >  0/  '  =  1/2,....  Then  from 

\ 

Theorem  2  one  has 

Corollary  5:  Assume  that  Z  satisfies  Condition  (Cl).  Then  VT  is  NBU, 
whenever  T.  is  NBU. 

i 

Note  that  when  the  components  are  independent  (Cl)  is  automatically 
satisfied.  Since  a  degenerate  random  variable  is  NBU,  one  can  easily  see 
that  Corollary  5  is  a  generalisation  of  Proposition  2  of  Ross  and  Schechtman 
(1979).  In  fact  the  above  result  holds  true  for  more  general  systems,  more 
general  than  those  considered  in  Chiang  and  Niu  (1980).  Such  a  generalisation 
is  discussed  in  Shaked  and  Shanthikumar  (1985). 


DEPENDENT  PARALLEL  SYSTEMS  WITH  PHASE  TYPE  SAFETY  PERIODS 


In  this  section  we  assume  that  (T^)*^  has  a  phase  type  distribution 
with  representation  (c,  C)  and  m  phases  (e.g.  Neuts  (1981)).  That  is  T.  has 
the  same  distribution  as  the  time  to  exit  from  a  lossy  process  (say) 

N  =  {N(t),  t  e  R+}  with  state  space  {1,2,...,m},  infinitesimal  generator  C 
and  initial  probability  vector  c.  Then 


and 


where 


E{T.}  =  c  C_11 


H(x)  =  c 


00  1  i  -nx  (nx) 

I  (I  +  n  C)  e  nx  - n -  1 

i=0  n 


r|  =  max  { -C( i , i ) }  . 

The  evolution  of  the  parallel  system  with  a  back-up  unit  that  has  a 
phase  type  lifetime  can  be  described  as  follows:  the  process  Z  evolves  as 
before  starting  from  some  initial  state  Z(0)  different  from  2-  Every  time 
Z  enters  state  \  an  auxiliary  process  N  is  initiated  into  state  i  with  probability 
Cj,  i  =  1,2,...,m.  The  process  N  is  lossy  and  is  governed  by  the 
infinitesimal  generator  C.  If  the  auxiliary  process  N  leaves  the  set  of  states 
{1,2,...,m}  before  Z  exits  state  2,  the  system  fails.  Otherwise,  the  auxiliary 
process  is  terminated  as  soon  as  Z  exits  state  2-  Let  Z  =  { Z ( t ) ,  t  e  R  +  }  be 
this  combined  process.  The  state  space  of  Z  is  then  S'U{(2,1  ),0,2),  •  •  •  ,0,m' 
Noe  that  a  state  0,0  means  that  Z  is  in  state  and  N  is  in  state  i,  i  =  1,2,. 

A  A 

The  infinitesimal  generator  Q  of  this  process  Z  is  given  by 


lx {1  > 2 , . . . , m j  _1‘.hCD 


-Ml  +  C 


An  example  with  Q  as  given  in  Figure  2  and  m  =  2  is  given  in  Figure  3. 
Note  that  Z  is  a  lossy  process  and  the  exit  time  is  V^..  Therefore 


E{Vt}  =  (a,  0)  Q'1  V  , 


P{VT  >  x}  =  (a,  0)|  1(1  +  1 


Q)'  e  px(0x)' 


■I1" 


x  >  0, 


where 


0  =  max  {~Q(i,i)} . 


Since  the  class  of  phase  type  distributions  is  weakly  dense,  one  may  use  the 
above  results  to  approximate  the  distribution  of  V-j.  when  the  safety  period 
do  not  have  a  phase  type  distribution.  We  will  next  consider  a  very  special 
case  of  this:  the  safety  period  is  a  fixed  constant. 


3.  SYSTEM  FAILURE  TIME  WITH  FIXED  SAFETY  PERIOD 


In  this  section  we  wll  assume  that  the  lifetime  of  the  back-up  unit  is  a 
fixed  constant  A.  That  is  T^  =  A  w.p.1,  j  =  1,2,....  Then  from  (5)  and 
(14)  one  obtains 


(0,0,0)  (1,0,0)  (0,1,0)  (0,0,1)  (1,1,0)  (1,0,1)  (0,1,1) 


and 


2i 


0 

E{ V A }  =  (a  +  eMA  b  -  b)  Q"1  V  +  eMA  JA  Mxe'Mxdx  +  A  . 
a  -  -  u  0 

The  probability  distribution  P{VA  i  x}  can  be  obtained  from  (7)  with 
appropriate  substitution.  Specifically, 

P{VA  i  x}  =  0,  x  ^  A 

and 

P{Vfl  S  x}  =  {eMA  I  (1-e'MA)j  F(j)  ®  G*(j)]  ®  Fn(x-A),  x  ^  A, 
A  j=0  U 

where  F  and  Fq  are  given  by  (14)  and 

G*(x )  =  (1-e‘Mx)/(1-e'MA),  0$x^A. 

In  the  spirit  of  Theorem  4  and  Equation  (12)  one  may  also  use  the  following 
exponential  approximations  for  P{VA  S  x}  when  e  MA  is  very  small. 

P{VA  <  x}  =  Fg  ®  B(x-A),  x  §  A  (17; 

and 

P{VA  ^  x}  S  1-exp{-(x-A)/(E{VA}-A)},  x  £  A, 

where 

B(x)  =  1-exp{-xe'MA/{(E{X)  +  E  {V*})  (1-e'MA)}},  (18 

and 

A 

E{Y*}  =  (1-e  )  J  M  x  e  dx. 

0 


In  (17)  we  approximate  the  time  to  first  failure  distribution  by  the 


'V"  \r*  r*v"  v  v  \-i 


convolution  of  a  phase  type  distribution  Fq  given  by  (14)  and  an  exponential 
distribution  given  by  (18).  Next  we  will  look  at  an  alternative  approximation 
for  P{VA  $  x}  . 

It  is  well  known  that  a  degenerate  distribution  can  be  arbitrarily  closely 

approximated  by  an  Erlang-m  distribution  with  sufficiently  large  m.  Let 

<T?r=i be  an  i.i.d.  sequence  of  safety  periods  with  Erlang-m  distribution: 

equivalently  phase  type  distribution  with  representation  (c  ,  C  )  where 

— rn  m 

cm  =  (1,0,0,. ..,0),  Cm(i,i+1)  =  m/A  =  -Cm(i,i),  i  =  1,2,...,m-1; 

Cm(m,m)  =  -  m/A  for  m  =  1,2, -  Note  that  E{T™}  =  A  for  all  m  =  1,2,... 

and  Var  {VI1}  =  A2/m,  m  =  1,2,....  Now  let  V-j-m  be  the  system  lifetime  VT 
obtained  from  (4)  with  ( T j ) j=-j  replaced  by  (Vj1)”.^.  Now  combining  the  above 
observations  with  the  results  in  Section  2  (Equation  (15)  and  (16))  one  obtains 


E{VTm]  =  (a,  0)  Qm'  V , 


P{V  m  £  x  }  =  (a,  0)j  I  (1+1  Q  ) 


ePrr  (Pmx)' 


i=0  ^ 


I1,  x  >  0, 


=  1,2,..., 


where 


=  r  Qo 

L  r  •  h(i 


-QoT-Sm  “ 

-Ml  +  Cm. 


and  8  -  max  {*Q  (i,i)}.  Thus,  for  large  m  one  has 


E  { V  A }  =  E{VTm} 


P{VA  i  x]  *  P{VTm  S  x}. 


*»  *» 


23 


REFERENCES 

Chiang,  D.T.  and  S.C.  Niu  (1980),  On  the  distribution  to  first  system  failure, 
Journal  of  Applied  Probability,  17,  481-489. 

Dunbar,  L.C.  (1984),  A  mathematical  expression  describing  the  failure  probability 
of  a  system  of  redundant  components  with  finite  maximum  repair  time, 
Reliability  Engineering,  7,  169-179. 

Keilson,  J.  (1979),  Markov  Chain  Models-Rarity  and  Exponentiality,  Springer- 
Verlag,  New  York. 

Keilson,  J.  and  W.R.  Nunn  (1979),  Laguerre  transformation  as  a  tool  for  the 
numerical  solution  of  integral  equations  of  convolution  type,  Applied 
Mathematics  and  Computation,  5,  313-359. 

Neuts,  M.F.  (1981),  Matrix-Geometric  Solutions  in  Stochastic  Models:  An 
Algorithmic  Approach,  The  Johns  Hopkins  University  Press,  Balitmore. 

Ross,  S.M.  (1984),  A  model  in  which  component  failure  rates  depend  on  the 
working  set,  Naval  Research  Logistics  Quarterly,  31,  297-300. 

Ross,  S.M.  and  J.  Schechtman  (1979),  On  the  first  time  a  seperately 

maintained  parallel  system  has  been  down  for  a  fixed  time,  Naval  Research 
Logistics  Quarterly,  26,  285-290. 

Schechner,  Z  (1984),  A  load-sharing  model:  the  linear  breakdown  rule, 

Naval  Research  Logistics  Quarterly,  31,  137-144. 

Shaked,  M.,  and  J.G.  Shanthikumar  (1985),  On  the  first  failure  time  of  a 
dependent  multicomponent  reliability  system,  manuscript  in  preparation. 

Shanthikumar,  J.G.  (1984),  Processes  with  new  better  than  used  first- 
passage  times,  Advances  in  Applied  Probability,  16,  667-686. 

Shanthikumar,  J.G.  (1985),  Bilateral  phase  type  distributions,  Naval  Research 
Logistics  Quarterly,  32,  119-136. 

Shanthikumar,  J.G.  and  U.  Sumita  (183),  General  shock  models  associated 
with  correlated  renewal  sequences,  Journal  of  Applied  Probability,  20, 
600-614. 

Sumita,  U.  (1981),  Development  of  the  Laguerre  Transform  Method  for 
Numerical  Exploration  of  Applied  Probability  Models,  Ph.D.  Thesis, 

University  of  Rochester,  Rochester. 


